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We study binary and the recently discovered process of ternary He-assisted recombination of 
ions with electrons in a low temperature afterglow plasma. The experiments are carried out 
over a broad range of pressures and temperatures of an afterglow plasma in a helium buffer gas. 
Binary and He-assisted ternary recombination are observed and the corresponding recombination 
rate coefRcients are extracted for temperatures from 77 K to 330 K. We describe the observed 
ternary recombination as a two-step mechanism: First, a rotationally-excited long-lived neutral 
molecule H3 is formed in electron-Hj]' collisions. Second, the H3 molecule collides with a helium 
atom that leads to the formation of a very long-lived Rydberg state with high orbital momentum. 
We present calculations of the lifetimes of H3 and of the ternary recombination rate coefficients for 
para and ortho-Hg". The calculations show a large difference between the ternary recombination rate 
coefRcients of ortho and para-H^ at temperatures below 300 K. The measured binary and ternary 
rate coefficients are in reasonable agreement with the calculated values. 



I. INTRODUCTION 

Electron scattering by the simplest polyatomic ion 
is of fundamental importance in plasma physics because 
H;^ ions are dominant in many hydrogen-containing plas- 
mas, including astrophysically-relevant plasmas in par- 
ticular. For quantum theory, the electron-Hjj' interaction 
is important because the process involves the simplest 
polyatomic ion that can be studied from first principles 
without adjustable parameters. In the long history of 
research on recombination of H;^ with an electron there 
have been numerous exciting, contradictory and some- 
times unaccepted results. We mention here only a few 
recent reviews devoted to recombination at low tem- 
peratures (see e.g. [U, 0, H, 13, Hi)- A successful theory 
of Hj^ recombination with electrons at scattering ener- 
gies below 1 cV was developed relatively recently when 
the Jahn- Teller non-Born-Oppenheimer coupling was in- 
cluded into theory 0, B Q • The measurements of the 
recombination rate constant in different storage ring ex- 
periments have also converged to the same value recently, 
after it became understood that the internal rotational 
and vibrational degrees of freedom of the Hj^ ion influence 
the recombination rate [1, [13, Ei • The theoretical devel- 
opments and improvements in storage ring experiments 
have resulted in a reconciliation of theory and experiment 
for the binary electron-H^ recombination process. 

However, the Hjj" and D^^ recombination experiments 
carried out in afterglow plasmas [1, 0, 0, 

M M M [H, 

[TgI . have repeatedly given rate coefficients very differ- 
ent from the ones obtained in the aforementioned storage 
ring experiments and the theoretical calculations. Until 
now, the plasma studies have not been fully understood 



and, as a result, they have been frequently rejected be- 
cause they do not mesh with the present understanding 
of the binary DR process (see, e.g., the very recent re- 
view discussing this subject [3). However, the experi- 
mental plasma results are reproducible and they demand 
an understanding and an integration into the full picture 
of H^ recombination with electrons. The present work 
discusses and, we hope, clarifies some aspects of this com- 
plex problem. 

The plasma measurements are usually carried out in 
a He/Ar/H2 gas mixture (see the review by Plasil et al. 
Q) or in a pure H2 gas [13, [11]. The main question is 
how to "reconcile" [l^| the rate constant (including its 
dependence on experimental conditions) observed in an 
H;^ dominated plasma with actual theory and with data 
from the storage ring experiments [l8| . The plasma ex- 
periments are typically carried out with helium buffer gas 
at pressures in the range 50-2000 Pa. It has been gen- 
erally accepted that such pressures are too small to pro- 
duce appreciable ternary helium-assisted recombination 
of Hj]'. The fact that the neutral-stabilized recombina- 
tion can sometimes play a role was predicted many years 
ago by Bates and Khare [2^ and confirmed for some ions 
(but not for Hj^) by experiments [U, The typical 
value for the three-body recombination rate coefficient 
A'hc with helium as an ambient gas obtained in exper- 
iment and estimated theoretically [2l|, is of order of 
10-27 cm^s-i at 300 K. Thus, at pressures of 1300 Pa, the 
corresponding apparent binary recombination rate coeffi- 
cient is smaller than lO-^cm'^s""'^, which would therefore 
be negligible in comparison with the now-acce pted bi- 
nary H3 recombination rate coefficient [1, ITol.llll. [lH . 
This would also be negligible in comparison with the bi- 
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nary rate coefficients (about 10~^ cm^s"^ at 300 K) for 
the majority of molecular ions. This is the reason why 
the ternary recombination of H;^ had previously been ne- 
glected in FALP (flowing afterglow) and SA (stationary 
afterglow) experiments carried out at 50-2000 Pa. 

In our recent study [1^ we have shown that at 260 K 
a significant fraction of the + c^ collisions leads 
to formation of long-lived (up to tens of picoseconds) 
rotationally-excited neutral Rydberg molecules Hg. The 
formation of long-lived Hg and collisions of Hg with he- 
lium atoms can influence the overall process of recombi- 
nation of the H^ dominated afterglow plasma. By mea- 
suring the helium pressure dependence of the recom- 
bination rate coefficient we have found that at 260 K the 
ternary recombination rate is comparable with the binary 
rate, already at pressures of several hundred Pa. The 
observed ternary recombination of ions is more effi- 
cient by a factor of 100 than the ternary recombination 
rate predicted by Bates and Khare This suggests 

that ternary recombination of H^ ions is associated with 
intrinsic features of the interaction between Hjj' and elec- 
trons at low collision energies p3| . The recombination 
process described by Bates [l^l is a ternary process of a 
different nature. 

This study presents a further extension of our measure- 
ments to a broader range of temperatures and it extracts 
the temperature dependence of the binary and ternary 
recombination rate coefficients of H^ ions. Before pre- 
senting these new results and their interpretation, some 
of us (the Prague group) wish to provide a few comments 
concerning the plasma experiments made in Prague in re- 
cent years. 

Previously, we have studied the H^ recombination in 
decaying plasmas formed from a He/Ar/H2 gas mix- 
ture using several afterglow experiments based on several 
modifications of flowing and stationary afterglow appara- 
tuses (FALP ^3, AIS A- Advanced Integrated Stationary 
Afterglow [l,lll,[l3|, TDT-CRDS-Test Discharge Tube 
with Cavity Ring-Down Spectroscopy [1^, [2^ ) . We have 
systematically investigated the dependence of the recom- 
bination process on hydrogen partial pressure. Initially 
(in 2000), our intention was to explain the influence of the 
internal excitation of H|^ ions formed in a sequence of ion- 
molecule reactions on the measured recombination rate 
coefficient. Experimental conditions were such that the 
plasma was decaying for a long time (up to 60 ms), so the 
internal excitation was quenched in collisions with He and 
H2 . As a rule, we skip the first 10 ms of the decay process 
considering it as a formation time. At low H2 densities, 
[H2] < 10^^ cm~^, and hchum pressures of 100-300 Pa, 
we have observed a decrease of the recombination rate co- 
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At high hydrogen densities, [H2] > 10^^ cm^"^, the mea- 
sured recombination rate coefficient was independent of 
[H2]. We will consider this again below when discuss 
equilibrium conditions for recombining H^ plasma. A 
similar dependence on D2 density was observed in a Dj]"- 



dominated afterglow plasma (see e.g. d, [i3)- 

In our early experimental publications [13, [131 1 "i*^" 
rived the recombination rate coefficient from experimen- 
tal measurements assuming that the process is strictly 
binary (R^ + e~) at the low-pressure limit. The as- 
sumption was based on the then-existing level of knowl- 
edge about recombination (see the recent and older 
reviews in Refs. i, S Q [H, [H, H^). However, we 
soon realized that the observed dependence on hydrogen 
density is coupled to the multistep character of the re- 
combination process in a plasma. The theory of binary 
dissociative recombination of is more immediately 
applicable to the storage ring experiments. Therefore we 
denoted the plasma recombination rate constant as an 
effective deionization rate constant aeff 0, [3 ■ Later ex- 
periments using absorption spectroscopy for ion density 
measurements and for the identification of the recombin- 
ing ions (TDT-CRDS) [H, [H, did not clarify the 
mechanism of recombination in an afterglow plasma. 



The interpretation of data from storage rin ; 
31 1 exp eriments and from afterglow 0. 
24 . HH 33[_experimcnts was at this point not reconciled 
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HjHiJsII- We stress here one principal difference be- 
tween the two types of experiments: In a storage ring 
experiment, the recombining ions and electrons have a 
small adjustable relative velocity; the measured cross- 
section corresponds to the binary electron-ion recombi- 
nation process. In a plasma experiment an ambient gas 
must be used. Typically, the pressure of the ambient gas 
in the afterglow experiments is 50-2000 Pa. The ions and 
electrons in afterglow plasma undergo multiple collisions 
with neutral particles (He and H2 in the experiments 
discussed here) prior to their mutual collisions; these are 
collisions at low energy (around ~ 0.025 eV). In a stor- 
age ring, collisions of Hj^ with neutral particles are also 
possible, but when an ion collides with a background gas 
molecule or atom, it is removed from the beam and does 
not contribute any further to the observable recombina- 
tion events Q . 

Plasma experiments measure the thermally averaged 
rate constant. It means that the afterglow plasma should 
ideally be in thermal equilibrium with respect to all de- 
grees of freedom, in particular with respect to the rota- 
tional, vibrational and nuclear spin states of Hj^. Molec- 
ular hydrogen present in the recombining plasma equi- 
librates the ortho/para-Hj^ ratio and should make the 
"kinetic" temperature of ions and electrons to coincide 
with the He temperature [13, [s^. An internal excita- 
tion of Hg ions obtained in the process of H+ for ma- 
tion [1, Q is also quenched in collisions with He atoms, 
but the population of the lowest energy levels (lowest or- 
tho and para-states) is governed by collisions with H2. 
One has to have in mind that proton transfer reactions 
(from ArH^ or HJ to H2) producing B.'^ occur at nearly 
the rate of elastic collisions, but the rate coefficients for 
reactions between H;^ and H2, which change ortho-H;^ 
to para-Hj and vice versa ("state changing collisions") 
arc a factor of 5-10 smaller, fcgc ~ 3 x 10"^" cm^s~^ 
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(see Refs. [3a, |33, H^). If recombination of is sensi- 
tive to the nuclear spin state of then recombination 
drives the ions out of ortho/para thermal equilib- 
rium. In contrast, the "state-changing collisions" with 
H2 shift the ortho/para distribution towards the ther- 
mal equilibrium. In this case, thermal equilibrium means 
a population of rotational states corresponding to ther- 
modynamic equilibrium at the given temperature. Note 
that we have already demonstrated that in a plasma at 
260 K the recombination of ortho and para-H^^ ions is 
very different psj . The recombination frequency (num- 
ber of events per unit time) is v^cc ^ Q^effi-e and the 
"state changing frequency" is Vsc ~ fcsc[Il2]. The condi- 
tion i^sc > i^icc required for thermal equilibrium leads to 
fcsc[Il2] > acsn-c- If values typical for FALP experiments 
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2 X 10^ cm"'^ and aeff ~ 1.5 x 10" 



are 



used in this inequality, we obtain the condition on the hy- 
drogen density: [II2] > 10^^ cm~'^. This is an important 
density value for the interpretation of data in the plasma 
experiments (see also the definition of the rate coefficient 
in plasma as it is discussed in [1^). It is presumably 
not accidental that the density 10^^ cm~^ is the same as 
the density at which the measured effective H;^ recom- 
bination rate coefficient changes its character. We have 
in mind that at temperatures below 300 K fcsc and acs 
are dependent on rotational excitation (ortho and para 
states) of II2 and than also the condition for thermal 
equilibrium is more complex. We will come to this point 
again later. 

Because a recombining plasma contains neutral 
molecules and atoms, collisions between ions and elec- 
trons are perturbed by the neutral particles, which can 
influence the observed effective plasma recombination 
rate. The perturbation could be significant if in a 
electron-ion collision a long-lived highly excited interme- 
diate neutral molecule is formed (see the discussion in 
Rcfs. @, d^). If the third particle is an atom from the 
buffer gas, the probability of such a collision is propor- 
tional to its pressure and, as a consequence, the apparent 
recombination rate coefficient will be pressure dependent. 
Similar phenomena are well known from studies of ion- 
molecule association reactions, which can have a binary 
(e.g. radiative association) and/or three-body character 
depending on the lifetime of the complex. The rate of 
these processes can depend on pressure and temperature 
(e.g. associative reactions, see the discussion in (40|). 
Third-body perturbation of binary recombination in 
a plasma was already mentioned in our recent study [2^ . 
In the present study, we show further experimental ev- 
idence for the phenomenon and further aspects of the 
theoretical interpretation, in conjunction with numerical 
calculation. 

The rest of the article is organized in as follows. First, 
we briefly describe the experiments in section |TT1 Then 
we present new experimental results in section Hill where 
we interpret the observed dependencies of the apparent 
binary recombination rate coefficients (aoff ) on hydrogen 
and helium pressure, and from that analysis, derive the 
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FIG. 1: Cryo-FALP. Plasma created in a microwave discharge 
is carried along the flow tube by helium carrier/buH'er gas 
(from left to right in the drawing). After addition of Ar (via 
the indicated gas inlet) the plasma is converted from He J 
dominated to Ar"*" dominated. Further downstream H2 is in- 
jected to an already-relaxed plasma (with thermalized elec- 
trons) and an H3" dominated plasma is formed. The electron 
density decay downstream from the hydrogen entry port is 
monitored by an axially movable Langmuir probe. 



binary and ternary recombination rate coefficients. Fi- 
nally, in section [TVl we introduce our theoretical descrip- 
tion of ternary recombination of an plasma, present 
our calculation of the lifetimes of the long-lived neutral 
molecules II3 formed in the electron-ion collisions, and es- 
timate the rate coefficient for the ternary channel of 
recombination. Section IVl summarizes our conclusions. 



II. EXPERIMENTS 



For measurements of pressure and temperature depen- 
dences of the recombination rate we have built the 
Cryo-FALP apparatus, designed to operate in the range 
77-300 K and at helium pressures adjustable from ~ 400 
to ^ 2000 Pa. UHV technology and high buffer gas pu- 
rity (the level of impurities is < 0.1 ppm) are used in the 
Cryo-FALP (Fig. [J). 

The Cryo-FALP apparatus is a low temperature high 
pressure variant of the standard FALP apparatus [41| . 
Here we will just describe the essential features of the 
new construction. In Cryo-FALP. plasma is created 
in a microwave discharge (10-30 W) in the upstream 
glass section of the fiow tube (at 300 K). Because of 
the high pressure, He'*' ions formed by electron impact 
then react in a three-body association reaction with 
two atoms of helium, and a HeJ dominated plasma is 
formed. Downstream from the discharge region, Ar is 
added to HeJ dominated afterglow plasma to remove he- 
lium metastable atoms (He™ in Fig. [5]) formed in the 
discharge and to form Ar"*" dominated plasma (see de- 
tails in 0, [13, [13 ) ■ Via the second entry port situated 
approximately 35 ms downstream from the Ar entry port, 
hydrogen (diluted in He) is introduced into the plasma, 
which at this point is already cold. Note that the po- 
sition in the flow tube is linked to the decay time. In 
a sequence of ion-molecule reactions, an Hj^ dominated 
plasma is formed shortly after the second entry port. 

We have used a numerical model to simulate the pro- 
cess of formation of Ht dominated plasma. Examples 
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FIG. 2: Left panel: Calculated dependences of ion for- 
mation and plasma decay in FALP after addition of Ar 
(7.7 X 10^^ cm~^) and H2 (5 x 10^^ cm"^). Argon was added 
in the flow tube 33 ms before H2. The time scale origin is 
set at the position of the hydrogen injection port. The right 
panel focuses on a narrow time interval corresponding to a 
transition from Ar^ dominated to Hj|" dominated plasma af- 
ter adding hydrogen. 



of ion density evolutions calculated within the model 
for conditions typical for the present Cryo-FALP exper- 
iments are shown in Fig. [2l Similar calculations were 
carried out for all of our experiments presented in this 
work. Argon also plays a role in plasma relaxation and 
in formation of ions. However, when an domi- 
nated plasma has already been formed, Ar does not play 
any role (in contrast with He) because its density is at 
least four orders of magnitude lower than the helium den- 
sity. Downstream from the Ar entry port the flow tube 
is thermally insulated and cooled to the desired temper- 
ature by liquid nitrogen. Because theory suggests a very 
strong temperature dependence of the process of inter- 
est, we monitored carefully the temperature of the flow 
tube. An axially movable Langmuir probe [i^ is used to 
measure the electron density decay downstream from the 
hydrogen entry port [2^ . Electron energy distribution 
functions (EEDF) were checked along the flow tube to 
characterize plasma relaxation during the afterglow time 
[H,!!!]. Recombination of ("benchmark ion" 
was used for calibration of the Langmuir probe over a 
broad range of pressures. 

An advanced analysis was used to fit the decay curves 
[1, [13, 13 with the purpose to obtain recombination rate 
coefficients from the measured decay of electron densi- 
ties. In the analysis we have assumed that once an 
dominated plasma is formed, the plasma decay can be 
described by a single value of the recombination rate co- 
efficient, which we call the effective apparent binary re- 
combination rate coefficient, aeff- The plasma decay can 
then be described by the balance equation with a recom- 
bination and a diffusion loss terms: 
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FIG. 3: Examples of electron density decay curves measured 
in the dominated plasma at several hydrogen densities. 
Upper panel: Cryo-FALP experiment. Lower panel: FALP 
experiment. The effective recombination rate coefficients ob- 
tained at different hydrogen densities are indicated. For com- 
parison, both panels also show the decay curves (rectangles) 
obtained in the Ar"'" dominated plasma with [H2] = 0. 



tively. We assume that plasma is quasineutral (no = n^)- 
The constant td characterizes the ambipolar diffusion 
during the afterglow. The recombination of H^ ions at 
temperatures below 300 K depends strongly on the ro- 
tational excitation of ions. The assumption about the 
constant value of aeff must be discussed for particular 
experimental conditions. 



In the present experiments we use helium densities in 
the range [He] ~ 6 x 10^^ - 6 x 10^'' cm~'^ and hydro- 
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where Uc and are electron and ion densities, respec- 



gen densities in the range [H2] ~ 10 
the Cryo-FALP and other experiments discussed we use 
normal hydrogen (n-H2) i.e. the mixture of ortho and 
para-H2 corresponding to 300 K (approximately 25% of 
para H2). The variation of the ortho/para-H2 ratio with 
temperature in the interval 77 K-300 K is not significant 
for the results of the present experiments. (In thermal 
equilibrium at 100 K the fraction of para H2 is ~ 38%, 
and at 77 K the fraction is ^ 50%). 
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III. EXPERIMENTAL RESULTS 

We monitored decay of the afterglow plasma in a 
He/Ar/H2 mixture at different temperatures and over 
a broad range of helium and hydrogen densities. Ex- 
amples of decay curves measured at 77 K and at 250 K 
at several hydrogen densities are plotted in Fig. [3l The 
dependence of the decay rate on hydrogen density is ev- 
ident. The obtained apparent recombination rate coeffi- 
cients (aoff) depend on all three parameters aoff(T, [H2], 
[He]). It clearly indicates that the observed "deionization 
process" is not pure binary dissociative recombination. 
In Fig. [3] we have also plotted the decay curves measured 
in a Hc/Ar afterglow plasma dominated by Ar"'' ions at 
77 K and at 250 K in otherwise very similar conditions. 
At 250 K (see lower panel) the decay curve is exponential 
(straight line in the semilogarithmic plot) because recom- 
bination of these atomic ions is very slow and the decay 
is governed by ambipolar diffusion. At 77 K (see upper 
panel of Fig. [3]) we observe a faster decay during the early 
afterglow (at higher electron densities). We assume that 
this faster decay is primarily due to coUisional radiative 
recombination (CRR) [1^, |4y] and partly also due to the 
formation of Ar J ions (in three-body association at low 
temperatures) followed by immediate recombination of 
these ions [43] . The rate of the decay is comparable with 
the rate calculated for the CRR at ~ 85 K. 

The apparent binary recombination rate coefficients 
(ttoff) obtained from measured decay curves in the H;^ 
dominated plasma are plotted as functions of hydrogen 
density in Fig. [H Examples of the data obtained in other 
experiments (FALP, AISA and TDT-CRDS) are also in- 
cluded in Fig. m Below we summarize the data plotted 
in Fig. m 

Upper panel: 

1. 77 K, Cryo-FALP. At a fixed flow tube temperature 
(T = 77 K) and at fixed [He] = 1.9 x 10^^ cm'^, 
the dependence of aes on H2 density was measured 
from [H2] ~ 2 X lO^^ cm'^ to ~ 10" cm'^. Then, 
for several other helium densities the recombina- 
tion rate coefficient was only measured in a lim- 
ited range of hydrogen densities close to [H2] ^ 
10^^ cm"-^. Wc plot only a few examples. 

2. 100 K, TDT-CRDS. The H+ ion density was 
measured by using laser absorption spectroscopy 
(CRDS technique). During the discharge pulse and 
during the recombination dominated afterglow, the 
ion temperature was determined from the Doppler 
broadening of an absorption line. The details of 
such experiments can be found in [ist . 

3. 330 K, TDT-CRDS. The absorption spectroscopy 
technique (CRDS) was used to measure the tem- 
perature and density of the recombining Hj ions 
[251 [26l . |48| . Relatively high hydrogen density was 
used in these experiments. Therefore, an extrap- 
olation had to be carried out in order to obtain 



tteff for lower hydrogen densities (see the discus- 
sion below). For details of this extrapolation see 
the discussion in [23, \^ . In some experiments, the 
He temperature was determined by measuring the 
Doppler broadening of the H2O line [i^. We have 
also assumed that the temperatures of H2O and 
He are identical. The H2O density is very low but 
sensitivity of CRDS is very high for this molecule. 

Lower panel: 

4. 130 and 230 K. AISA. Examples of data measured 
with AISA. The data were measured at a fixed tem- 
perature and at a fixed pressure as a function of 
hydrogen density. 

5. 170, 195, 250 K, FALP. Three different versions of 
FALP designed to work at high He pressures were 
used in these experiments. 

6. 300 K, Experiment by Laube et al. [i^. Low 
pressure FALP experiment: [H2] = 2x10^** cm""^, 
[He] = 1.6 X 10^^ cm~^, obtained aoff is 7.8 x 
10"® cm-'^s^^ 

7. 295 K, experiment by Gougousi et al. [l3l. In this 
FALP experiment the dependence of the recombi- 
nation rate coefficient on [H2] was observed. The 
helium pressure was about 130 Pa. 

8. 300 K, theory The value calculated for bi- 
nary dissociative recombination (for [H2] = 0). 

Notice that in certain experiments the rate coeffi- 
cients were measured over a limited range of hydro- 
gen densities. In our experiments, we have covered the 
101° cm-3 < [H2] < 10^6 cm-3 range (AISA, FALP, 
Cryo-FALP, TDT-CRDS). 

In the measured dependences shown in Fig.[3|wc distin- 
guish three clearly different regions of hydrogen densities 
that exhibit specific behavior of aoff as a function of [H2] . 
We indicate these regions as < 1, iV > 1, and N ^ 1. 
The "vertical shifts" for some dependences in Fig. [3| (such 
as the data presented with open and full triangles in the 
lower panel) are discussed below. Wc characterize the 
three regions as follows. 

1. [H2] < 10^^ cm^^, < 1. At these hydrogen den- 
sities ttoff decreases with decreasing hydrogen den- 
sity. At such conditions the Hj^ ions formed by 
an exothermic proton transfer (from ArH^ or Hj 
to H2) do not have enough collisions with H2 to 
establish ortho/para-H^ thermal equilibrium prior 
to their recombination (see e.g. the discussion in 
m, [13). The number of these H2 + H+ 
collisions that a Hj^ ion will undergo prior to its 
recombination with an electron (at typical condi- 
tions of the discussed afterglow experiments) is in- 
dicated in the Fig. [H by number N. li N < 1, 
the decay of the plasma caused by the recombina- 
tion (at a given electron density) is faster than the 
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FIG. 4: Measured dependencies of the effective recombina- 
tion rate coefficient of ions on fiydrogen density. The 
data obtained by Cryo-FALP and TDT-CRDS are plotted in 
the upper panel. The data obtained in other experiments are 
plotted in the lower panel. We have also plotted examples 
of data obtained in our previous studies using: AISA, TDT- 
CRDS and FALP. The FALP data by Laube et al. ^ and 
Gougousi et al. [T^ are plotted for comparison. Additionally, 
we show the value of the theoretical recombination rate coef- 
ficient calculated for binary dissociative recombination (DR), 
Q:Biii(300 K) [l,[3l|. Note that for this theoretical result, only 
the binary DR process is considered, and as such it cannot 
depend on the hydrogen density. 



rate of rcthcrmalization. Therefore, in this non- 
equihbrium regime, the individual state composi- 
tion of H;^ afterglow plasma is different along the 
flow tube because the effective recombination rate 
depends on the absolute value of electron density, 
which varies along the length of the tube. By "state 
composition" we mean not only the kinetic energy 
distribution, which is established in collisions with 
He atoms with a nearly collisional rate, but also 
the rotational and ortho-H;^/para-H J state distri- 
bution. A quantitative description of this particu- 
lar regime would require a much deeper theoretical 
analysis. 



2. 10^2 < [H2] < 5 X 10^3 > 1. 

In this regime, the measured rate coefficients are 



nearly independent of [H2]. On the basis of the 
same arguments mentioned above, it is clear that 
at [H2] > 10^^ cm^'^, the H;^ ion (formed shortly af- 
ter hydrogen is injected) has several collisions with 
H2 prior to its recombination with an electron. Be- 
cause only some collisions lead to a change in ro- 
tational excitation of H;^ or in ortho<->para transi- 
tions ^MilSZ]): we assume that if [H2] > 10^^ cm~'^ 
the recombining ions in the flow tube will be in 
thermal equilibrium corresponding to the hydrogen 
temperature, which is assumed to be equal to the 
helium temperature. Because of the independence 
of aoff on [H2] we will call this region the "sat- 
urated region". The boundaries of the saturated 
region depend on actual helium pressure, temper- 
ature and electron density. The region is broad 
enough to find conditions where the value of Ofcff is 
nearly constant (the plateau part of the Q;cfr([H2]) 
dependence). Depending on experimental condi- 
tions we have covered either the whole saturated 
region or in some cases just a part of it. 

3. [H2] > 5 X 10^^ cm"^, TV > 1. Here the mea- 
sured recombination rate coefficient increases with 
hydrogen pressure. This behavior is caused by a 
formation of H^ ions and their fast recombination. 
The process is temperature and pressure depen- 
dent (see details in [H, \E^, [5l|, [12]). Because of 
the strong temperature and pressure dependence 
of ternary association, the onset of this region de- 
pends on these parameters. 

At first sight the "vertical shifts" of the dependences 
plotted in Fig.[4]are very chaotic. We demonstrate below 
that they arise because the apparent binary recombina- 
tion rate coefficient (aoff) depends not only on hydrogen 
density and temperature but also on the helium density. 
In addition, we also show that the temperature depen- 
dence of Qicff is not monotonic. 

We will not discuss the < 1 region here. How- 
ever, in connection with the data obtained at 77 K we 
want to point out one difference from the data obtained 
at higher temperatures. In all high temperature experi- 
ments made at [H2] < 10^^ cm^'^, we have observed a fast 
decrease of aoff with decreasing [H2] (see e.g. the FALP 
data measured at 250 K in the lower panel of Fig. [3]). 
In measurements at 77 K using Cryo-FALP (see the up- 
per panel of Fig. [4]) the decrease is substantially smaller. 
The difference can be partly associated with the effect 
of collisional radiative recombination (CRR) at 77 K. At 
[H2] ~ 10^^ cm~'^ when the overall recombination rate 
coefficient aoff (77 K) > 1.0 x 10^^ cm'^s^^ we can neglect 
the influence of these processes (see the upper panels of 
Fig. Eland Fig. HI). 

Figure 3] has a great deal of information, because it ac- 
tually shows aoff as a function of three variables T, [H2], 
and [He]. By analyzing the data we found the form of 
the function aoff (T,[H2],[He]) in the "saturated region", 
i.e. for the plasma where the H'^ ions before recombining 
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surement of the pressure dependence at a fixed tem- 
perature. The straight line is the best fit through 
the measured data. 

2. 100 K, TDT-CRDS. The values plotted were ob- 
tained from the measured dependences of Ofefi on 
hydrogen density as limits approaching the satu- 
rated region (see upper panel of Fig. |4]). More de- 
tails are given in [53|. 

3. 100 K, Cryo-FALP. The present measurements. 

4. 130 K, AISA. The data were obtained from the 
measured dependence of acs on hydrogen density 
(see lower panel of Fig. |4]) . More details is given in 
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FIG. 5; The effective binary recombination rate coefficients 
(aeff) measured at the stated temperatures are shown as func- 
tions of the helium density. Upper panel: Low temperature 
data - Cryo-FALP (77 K), TDT-CRDS & Cryo-FALP (100 K) 
and AISA & FALP (130 K). The horizontal line indicates the 
theoretical value for binary dissociative recombination at 77 K 
[1,|38(. For details, see the summary in the text. Lower panel: 
Higher temperature data - Cryo-FALP (305 K), FALP (250- 
260 K, 170 and 195 K), TDT-CRDS (330 K), AISA (250-260, 
230 and 150 K). Individual points measured in other laborato- 
ries: Smith and Spanel [l^, Laube et al. [i^, Leu et al. [sil] 
(see description in the text). The horizontal line indicates the 
theoretical value for the binary dissociative recombination of 
H+ at 300 K [i,[33. 



undergo a sufficient number of "state changing collisions" 
with H2 to reach thermal equilibrium. For a better anal- 
ysis of the experimental data we plotted aoff measured 
at a fixed temperature as a function of helium density. 
The rate coefficients measured in the three principally 
different afterglow experiments at several different tem- 
peratures at hydrogen densities corresponding to the sat- 
urated region arc plotted as functions of helium density 
in Fig. El 

We briefly summarize the data plotted in Fig. [5] Upper 
panel: 

1. 77 K, Cryo-FALP. The data from measurements at 
[H2] ~ 10^^ cm~^ at different [He] (see upper panel 
of Fig.[?]). The experiment was intended as a mca- 



5. 130 K, Cryo-FALP. The present measurements. 

6. 77 K Theory. The horizontal line indicates the the- 
oretical value of the recombination rate coefficient 
for binary dissociative recombination (DR) of H;^ 
at 77 K 

Lower panel of Fig. [5) 

7. 150, 230 and 250 K, AISA. The values shown were 
obtained from measured dependencies of a^s on hy- 
drogen density as limits approaching the saturated 
region (see upper panel of Fig. H] and Ref. [H]). 
The accuracy of the values is high because the val- 
ues were obtained from a number of measurements. 

8. 170 and 195 K, FALP & AISA. The FALP data 
were obtained similarly to the AISA data, i.e. 
as limits approaching the saturated region. The 
straight lines connect the measured FALP points 
with the AISA points, which are obtained by an ex- 
trapolation of the data measured by AISA at 130- 
230 K. 

9. 250-260 K, AISA & FALP. Compilation of data 
from several experiments (for details see Refs. [2^, 

Hi). 

10. 305 K, Cryo-FALP. The data were obtained by 
directly changing the helium pressure in the flow 
tube. 

11. 330 K, TDT-CRDS. The values were obtained as 
limits approaching the saturated region (see upper 
panel of Fig.H]). 

12. 210 and 300 K, Leu et al. [5l|. In the experiment a 
microwave technique was used to measure the elec- 
tron density. The values 2.3 x 10~^ cm'^s"^ and 
3.3 X 10"'^ cm^s-i for 300 K and 210 K respec- 
tively (pressure ~ 2.6 kPa) were taken from Figs. 
2 and 4 of Ref. [U. 

13. 300 K, Laube et al. & FALP (FALP-MS) ex- 
periment. The measured value is aoff = 7.8 x 
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10~® cni^s"^ (see Fig.|4|). He pressure was ^ 70 Pa. 
Used hydrogen density corresponds to the satu- 
rated region. 

14. 300 K, Smith & Spanel ^ a FALP experiment. 
We show the value from their plot of the recombina- 
tion rate coefficient 0jS &- function of position along 
the flow tube (in Fig. 01). In the conditions that 
arise shortly after injection of hydrogen, Smith & 
Spanel measured the value aeS ~ 6 x 10~^ cm^s~^ 
for a relatively long time (see Fig. [4] in [IB])- Fur- 
ther downstream, they obtained a lower value of the 
recombination rate coefHcient. The helium pressure 
was ^ 260 Pa. 

15. 300 K, Gougousi et al. [ll]. We did not include 
their value in the graph because their measure- 
ments were at hydrogen densities too high and 
therefore out of the saturated region. We only men- 
tion their experiment in order to show that we in- 
cluded this study in our considerations. 

With regard to data obtained by Laube et al. [iot . 
Smith & Spanel [ll|, and Gougousi et al. [l3|, we have 
not analyzed their experiments in full detail. However, 
it is clear from their studies that, in agreement with our 
experimental data, there is a general trend: The effective 
rate constant a^fi increases with the increase of helium 
density (in the temperature interval covered). In Refs. 
[isl . [igI . [4^ the authors used relatively low helium pres- 
sures. Therefore, the increase of aofi with [He] was not as 
large as we observe in our measurements, such as Cryo- 
FALP. Very high pressure was used in [51|. As a result, 
they obtained large acfi in agreement with the trend. 

The experimental data plotted in Fig. [5] show that the 
apparent effective binary recombination rate coefRcient 
aeff measured at a fixed temperature depends linearly 
on helium density. We will discuss a possible recombina- 
tion mechanism below, but at this point we can assume 
that the process has a binary character at very low [He] , 
whereas with increasing [He] the helium assisted ternary 
process begins to contribute substantially to the overall 
recombination deionization process. Therefore, we can 
write for the observed linear dependence 



aoff = aB\n{T) + /VHc(r)[He] 



(2) 



in terms of the rate coefficient aein {T) for binary recom- 
bination and the rate coefficient K^^dT) for ternary He- 
assisted recombination. The two coefficients depend on 
temperature. Note that for some temperatures we have 
the FALP data (e.g. 77 K), for others we used AISA & 
FALP data, and also the TDT-CRDS data. The AISA 
data were obtained at low helium densities. The FALP 
data were taken at high helium densities and in most 
cases we could vary the He pressure in the FALP exper- 
iments. We extrapolate the data from AISA to temper- 
atures 170 K and 195 K. Then we plot a straight line 
through these new points (obtained by the extrapola- 
tion) and the corresponding points measured by FALP 




7-(K) 



FIG. 6: The binary recombination rate coefficient mea- 
sured in the present study (circles). The theoretical values 
shown (the dashed line) are calculated for relative populations 
of para and ortho corresponding to the thermal distribu- 
tion at the stated T 0, [s^ . The dotted line indicates values 
calculated from cross sections obtained in storage ring exper- 
iments UH, using cold ion sources. The present values 
are obtained as limits approaching the saturated region at low 
helium density from plots in Fig. [5] We also show the data 
(squares) obtained by Amano in pure hydrogen [ij. Further 
details are given in the text. 



(see lower panel in Fig.O. Using the straight line fits we 
obtain experimental values for aein and A'ho at 170 K 
and 195 K. 

The obtained values of the binary recombination rate 
coefficients aBin(r) are plotted in Fig. [6] as a function of 
temperature. We also show on Fig. [6] thermal rate coeffi- 
cients calculated for binary dissociative recombination of 
[1, HI]- The agreement is very good. Note that the 
data for 170 K and 195 K (at high [He]) were published 
Hi, mi before the calculation [|. The data for 250 K 
was partially obtained 3 also before the calculations. 
The pressure dependence of the 250 K data set was mea- 
sured after [2^ the calculations were published. Some 
data for 100 K and 300 K were also measured earlier, 
using CRDS 

Figure [S] also presents the data obtained by Amano 
using absorption spectroscopy [l3 | in experiments made 



with pure hydrogen at ~ 40 Pa ([H2]'^ 10 cm~^) used 
as a buffer gas. The large values of the measured re- 
combination rate coefficients indicate that the molecular 
hydrogen is a more effective three body partner (with an 
effective ternary rate constant of order of 10"^'^ cm^s^^) 
than helium, which is not surprising because H2 has ro- 
tational and vibrational degrees of freedom. 

We have also checked that the observed linear pres- 
sure dependence of ttes is not associated with the Lang- 
muir probe operating at different pressures from 40 to 
2600 Pa. For this test we have used the same probe to 
measure the recombination rate coefficient for ions 
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FIG. 7: The measured ternary recombination rate coefficient, 
KHe{T), for helium assisted recombination of Hg" with elec- 
trons. The legend indicates the experiments used to extract 
the data. The plotted line is shown merely to guide the eye. 



[3, 01 • Our assumption was that q;q+ is pressure in- 
dependent because recombination is a direct binary 
(electron-ion) process. We have also studied the pressure 
dependence of recombination rate coefficients for HCO^ 
and DCO"'' ions [s^l using the probe. For both cases 
(Oj and HCO+ZDCO"*") we observed no pressure de- 
pendence. Our HCO^ and DCO"'' rate coefficients are 
in good agreement with results of Leu et al. [5l[ and 
Amano [14| measured by different techniques. Leu et al. 
[Slj used a microwave technique to determine electron 
densities but obtained results are consistent with our ob- 
servation. 

The ternary rate coefficients obtained from the data 
plotted in the Fig. O arc presented in Fig. [7] as a func- 
tion of buffer gas temperature. The values obtained for 
300 K are: aBm(300 K) = (4.7 ± 1.5) x IQ-s cmh-\ 
ifHolSOO K) = (2.5 ± 1.2) X 10-25 cm^s-i. The figure 
shows that the ternary rate coefficient has a maximum at 
~ 170 K. Towards lower temperatures the rate coefficient 
is decreasing. This is a surprising result if one takes into 
account studies of ternary association processes. For such 
processes the ternary rate coefficients decrease monoton- 
ically with temperature (see Refs. [sHIstI. [sst ) . 

We briefly summarize the data plotted in Figure [71 



1. 77 K, Cryo-FALP. The data obtained from the pres- 
sure dependence measured in the present experi- 
ments at 77 K. See upper panels of Figs [31 [H and 
[51 The value of i^He was obtained from measure- 
ments that were repeated many times. 

2. 100 K, Cryo-FALP & TDT-CRDS. The data ob- 
tained by a combination of the values measured 
in two experiments. The TDT-CRDS values are 
based on measurements of the ion density evolu- 
tion during the afterglow using CRDS absorption 



spectroscopy, (see the upper panels of Figs. [Hand 
[5]). In order to calculate the ion density from the 
absorption signal, thermodynamic equilibrium was 
assumed. The assumption should be valid at the 
hydrogen densities used. The Cryo-FALP values 
were measured in the present experiments (upper 
panel of Fig. [5]). 

3. 130, 170, and 195 K, AISA & FALP. Data obtained 
in two experiments (Fig. [5]). 

4. 150 and 230 K, AISA. AISA was used to measure 
the dependence of acs on hydrogen density over 
a broad range including the "saturated region", 
(lower panel of Fig. [J). From these measurements, 
average values of the recombination rate coefficient 
in the "saturated region" have been obtained for 
several temperatures (see Fig. [5]). These data are 
very accurate because the measurements have been 
done many times at a fixed temperature. As we 
did not measure the dependence on helium pres- 
sure in the AISA experiments, we calculated 
using Eq. ^ with the current theoretical value 
for Q;Bin(7')- The data have been obtained at low 
pressures, 160-330 Pa. Therefore, the ternary rate 
coefficients extracted from this data have large er- 
ror bars. 

5. 250-260 K, FALP & AISA. These data were ob- 
tained in several experiments. A detailed descrip- 
tion was given in our previous publication |23| but 
there is one key difference: In that study, [2^ the 
"260 K" plot also included values obtained at tem- 
peratures close to but different from 260 K; the 
data shown there for Q!off(260 K) were recalculated 
from the measured values, assuming a T~^'^ depen- 
dence, which is a small correction. Because of the 
additional experiments with Cryo-FALP the recal- 
culation procedure was not necessary in the present 
work. Now it is clear that acff depends on both 
temperature and pressure; in the vicinity of 260 K 
the temperature dependence is steeper than T~'^'^. 

6. 305 K, Cryo-FALP. In this experiment the pressure 
dependence was measured (see Fig. [S|). 

7. 330 K, TDT-CRDS. The measurements were sim- 
ilar to those at 100 K (upper panel of Fig. [H and 
lower panel of Fig. [5]). Two different absorption 
lines were used in these studies. We have obtained 
the same value of the rate coefficient a^ff using ei- 
ther absorption line. 

8. 210 and 300 K, Data from Leu et al (see 
Fig. [5]). We have utilized the current theoreti- 
cal binary recombination rate coefficients for 210 
and 300 K and Eq. [21 to obtain the corresponding 
ternary rate coefficients. 
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IV. THEORETICAL MODEL FOR 
HELIUM-ASSISTED NEUTRALIZATION OF 
THE AFTERGLOW PLASMA 

It is possible to estimate theoretically the rate coef- 
ficient, Eq. (21), of the He-assisted recombination of 
with electrons in the following way. To stress that it is an 
estimated value, we use symbol K^'^ for the theoretical 
coefficient. 

He-assisted recombination is treated as a two step pro- 
cess. In electron-Hjj" scattering, rotational autoionization 
resonances of the neutral molecule Hg play an important 
role (see, for example Ref. 0, HI])- The star next to H3 
refers to the unstable character of these autoionizing res- 
onances. Such resonances with angular momentum I ~ 1 
could be very quite broad due to the high probability 
to capture an electron into the I = 1 partial wave, but 
they normally contribute relatively little to the two-body, 
H;^ -|- e~ recombination 0|, because almost always, such 
resonances decay back into a free electron and an H;^ 
ion. As we demonstrate below, lifetimes of some of these 
resonances can be quite long. If, during their lifetime, 
the Hg molecule collides with a helium atom, the colli- 
sion can lead to a change in the electronic state of the 
outer electron or in the rotational state of the H^^ ion and 
therefore, can make the autoionization process impossi- 
ble (or, at least, much less probable than the dissociation 
of H3). The overall rate coefficient for such He-assisted 
recombination is given by the formula derived in Ref. 



if 3d = k,At( 



(3) 



3 ' 



where a* is the rate coefficient for the formation of H 
ki is the rate coefficient for Z-changing collisions between 
He and Hg leading to the eventual dissociation of Hg. 
Finally, At is the delay time in the Hj^ -I- c~ collision in 
the sense introduced by Smith [G^I): it is an additional 
time that the electron spends near the ion in the modi- 
fied Coulomb potential compared to the coUisional time 
in a pure Coulomb potential. The delay time and the 
coefficient a* are substantially different from zero only 
near resonances. Therefore, the three-body rate coeffi- 
cient discussed above is expected to vary resonantly as a 
function of coUisional energy. 

We stress here that the rate constants ki, and a, de- 
scribed above depend on the corresponding scattering en- 
ergies (and defined as usual as cross-sections multiplied 
by relative velocities). They are not yet thermally av- 
eraged over the Maxwell-Boltzmann distribution. Cor- 
respondingly, the ternary rate constant K^'^ in Eq. ([3]) 
depends on two energies: the energy E of the + e~ 
collision due to the dependence of At and a, on E; and 
the energy Eb_c of collision between He and Hg due to the 
dependence of ki on E'hc- In this approach, however, the 
coefficient fc;(i?Ho) is considered to be constant over the 
energy interval of interest. The rate coefficient a* in our 
estimation is evaluated as va{E), where v ~ (2£'/m)^/^ 
is the relative velocity, <t{E) is the cross-section for the 



process leading to the delay time, m is the reduced mass 
of the + e~ system. 

For the following discussion, we assume that at a given 
energy E, there could be several open H;^ -I- e~ ioniza- 
tion channels [i = 1,2... rio). Such a situation is possi- 
ble for the conditions of the present experiment. If there 
are several open channels, the three-body coefficient K^'^ 
above should be averaged over the incident channels and 
summed over the final ones. If the incident ionization 
channel (before collision) is denoted by the index i, and 
the final one (following the collision) is j, then the corre- 
sponding rate coefficient for the three-body recombina- 
tion during the i — > j inelastic collision is given by 



kiAtj^y/2E/maj,{E) 



(4) 



with <Jji{E) being the cross-section for the inelastic col- 
lision. The delay time Atji for the process is an element 
of the delay-time matrix At as it is introduced and dis- 
cussed by Smith [60l |: 



At-, 



SJ^ dE 



(5) 



where Sji is an element of the scattering matrix for the 
i j process. Notice a difference in conventions in the 
present paper and Ref. (60j : Here, the first index in 
each matrix denotes the final channel. The second index 
denotes the incident channel. In Ref. [G^I , the convention 
adopted was the opposite, i.e., the first index ~ incident 
channel, second index ~ final one. 
The cross-section aji{E) is given by 



a,,{E) = 



(6) 



Using equations (jH El O and taking the sum over the 
final ionization channels j, we obtain the three-body rate 
coefficient Kf^ if the initial state of the + e~ system 
is i: 
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The sum in the second line can be simplified as [6C 
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where Qa is a diagonal element of the lifetime matrix 
introduced by Smith (60j 
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(8) 
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FIG. 8: Diagonal elements Qa of matrix Q for the three low- 
est (rotational) incident channels for the e~ + collisions. 
The rotational channels are {N+,K+) = (11), (10), and (22). 
Each maximum in Qa corresponds to an autoionization reso- 
nance. The lifetime of a resonances is given by Qu/i eval- 
uated at the maximum if there is only one channel open, 



and the product S'^ ^ in the above equation is the regu- 
lar matrix product. Due to the aforementioned difference 
in matrix index conventions, the order of the product is 
opposite to the one in Ref. Therefore, for Kf'^ we 

have 



3d 



or in atomic units 



(9) 



(10) 



The next step in the evaluation of the three-body rate 
coefficient is the thermal average over incident ionization 
channels i and over the Maxwell-Boltzmann velocity dis- 
tribution for a given temperature T. The average over 
incident channels is given by 



E,w,exp(-^ 



(11) 



where Wi = {2N^ + l)(2/i -I- 1) is the statistical weight of 
the incident channel, A^^^ and Jj are the angular momen- 
tum and the nuclear spin of the Hj^ ion respectively, Ei 
is the energy of the incident channel (rotational energy 
of channel i) . The average over the velocity (energy) dis- 
tribution should be performed over the two energy vari- 
ables, E and i?Ho, in Kf'^iE, Ehc)- Because Kf'^ does 
not depend on Enc the thermal averaging is reduced to 
the familiar two-body averaging integral over E only 



^7r(fcBT)3 Jo 



/vf (£;)exp I - 



E 



EdE . (12) 



Combining the above equations, we obtain 

2 ir Kf'w.^ exp (-^) VEAE 



yf^^f^BTfY.^ exp (-tIt) 
or using the lifetime matrix element Qa^ 



2tt ^iT^ilo QiiWieyiY> 
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(14) 
from 
In 



The lifetime matrix Q is calculated using Eq, 
the scattering matrices obtained numerically 
practice, the scattering matrix S^^-* (and matrix Q^^^) 
are calculated for a fixed value of the total angular mo- 
mentum N of the Hjj" -|- c~ system, which is the sum 
of the ionic TV"*" and electronic I angular momenta. The 
contributions Q\f^ from different A'' should be accounted 
in Eq. (fH|) . namely as 



N 



(AT) 



(15) 



The sum in the above equation runs over all N for which 
the incident channel i enters into the scattering matrix. 
Since the principal contribution to the rotational captur- 
ing of the electron comes from the I = 1 partial wave, the 
sum has three or less terms. 

The coUisional /-changing process is known to be rel- 
atively effective for excited atomic Rydberg states [6l|. 
Because the = 1) molecules formed in the H;^ -I- e~ 
collisions have a large principal quantum number (n ^ 
40-100, see Fig. [5]), it is reasonable to use the atomic 
rates for Z-changing collisions in our estimates. For 
Na* -f- He /-changing collisions the experimental rate 
is 2.3 X 10^^ cm'^s^^. In our estimation, we use this 
value for ki. Using the procedure described above and 
the value above for ki we have calculated thermally av- 
eraged rate coefficient for the He-assisted recombination 
of . The coefficient is shown in Fig. [5] as a function 
of temperature. The overall agreement with the experi- 
mental rate coefficient shown in Fig. [7] is reasonable given 
the approximative approach that we used here. Below 
150 K, the experimental curve drops down, but the the- 
oretical (dashed) curve continues to grow. Interestingly, 
the curve for pure ortho-HJ^ (dot-dashed curve) has a be- 
havior similar to the experiment. In fact, it is possible 
that in the experiment the ortho and para-H^ are not 
in thermal equilibrium at low temperatures. Our previ- 
ous calculation Q showed that the binary recombination 
rate of H;^ with electrons is much slower for ortho-H;^ 
than for para-H;^: At low temperatures, due to fast re- 
combination of para-H;^ the ratio ortho/para-H;^ could 
be larger than the one at thermal ortho/para equilib- 
rium. In such a situation the averaged theoretical rate 
coefficient K^'^ in Fig. [5] would be closer to the curve for 
ortho-H|^, i.e. it would be lower at small T similar to the 
experimental behavior (Fig. [7]). 
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FIG. 9: Calculated thermally-averaged three-body rate co- 
efficient (^K^'^y The rate coefficients calculated separately 
for ortho and para-H^ are very different. If the recombin- 
ing plasma is not in thermal equilibrium with respect to or- 
tho/para ratio, the averaged rate coefficient (dashed) could 
be very different from the one shown. 

V. CONCLUSIONS AND DISCUSSION 

We have studied recombination of ions in an after- 
glow plasma, in the presence of a helium buffer gas with 
a small admixture of molecular hydrogen. The helium 
densities were in the range 0.5 - 6 x 10^^ cm"'^ and hy- 
drogen densities 1-100 x 10""^^ cm""^. In such conditions 
the H;^ ions formed in the plasma have several collisions 
with H2 before they recombine with electrons. Thus, 
we assume that in these conditions the ions are in or- 
tho/para thermal equilibrium. The apparent binary re- 
combination rate coefficient a^s was measured as func- 
tion of hydrogen and helium densities for several tem- 
peratures in the 77-330 K range. From the experimental 
data we have derived the binary and ternary recombina- 
tion rate coefficients and their dependences on tempera- 
ture. The measured binary recombination rate coefficient 
is in good agreement with recent theoretical calculation 
over the whole covered temperature range (77-330 K). 
Therefore, for the first time, the recombination rate co- 
efficients obtained in afterglow plasma experiments agree 
with storage ring experiments and with theoretical val- 
ues. As we have demonstrated in the present study, re- 
sults from previous afterglow plasma experiments were 
previously interpreted without taking into account the 
role of the buffer gas and, as a result, those experiments 
seemed to disagree with the storage ring experiments and 
with theoretical calculation. The present work reconcile 
observation data from the plasma and storage ring ex- 
periments and the theoretical result. 



The obtained binary rate coefficient at 300 K is aBin = 
(4.7±1.5) X 10~* cm^s~^. The observed ternary recombi- 
nation (iCHe(300 K) = (2.5± 1.2) x lO^^s cni^^s-i) is fast 
and at pressures of hunderds of Pa it is already dom- 
inant over the binary process. The dependence of the 
measured ternary recombination rate coefficient on tem- 
perature has a maximum at ~ 130-170 K. The observed 
ternary process is more effective by factor about 100 than 
the process of ternary recombination described by Bates 
and Khare 

To explain the process of fast ternary recombination 
we have developed a theoretical model for the process. In 
particular, we have calculated the delay time in + e~ 
collisions (At as introduced by Smith ^^). We found 
that the delay time is sensitive to the rotational and 
nuclear-spin states of the H;^ ion. The delay time at 
collision energies E ^ 150 cm^^ can be of order of 
100 ps for para-Hj^. During that collision time, the H3 
molecule (H^ + e~ complex) can collide with a helium 
atom, which enhances the overall plasma neutralization. 
The calculated delay time was used to derive the ternary 
recombination rate coefficient. The derived ternary coef- 
ficient as a function of temperature is smaller than the ex- 
perimental value by a factor of order 2-10, which is plau- 
sible agreement for such a sophisticated process (from 
a theoretical ab initio point of view) within the some- 
what heuristic theoretical method employed. Theory can 
in principle be further improved. In particular, the co- 
efficient for Z-changing collisions can be evaluated more 
accurately. At temperatures below 130 K there is a qual- 
itative difference between measured and calculated val- 
ues of the ternary rate coefficient: the experimental rate 
constant decreases with temperature, the theoretical one 
continues to grow. It is possible that at low temperature 
ortho and para-H;^ ion are not in thermal equilibrium in 
this afterglow plasma because of the very different binary 
rate constants aBin(T) that have been predicted for low 
temperatures. Nevertheless, for a first semiquantitative 
study of this kind, the agreement between theory and ex- 
periment for the ternary rate coefficients is encouraging. 
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